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Abstract 

The concept of energy- dependent forces in quantum mechanics is re-analysed. 
We suggest a simplification of their study via the representation of each self- 
adjoint and energy-dependent Hamiltonian H = H{E) with real spectrum 
by an auxiliary non-Hermitian operator K which remains energy-independent. 
Practical merits of such an approach to the Schrodinger equations with energy- 
dependent potentials are illustrated using their quasi-exactly solvable example. 
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1 Introduction 



An introduction of phenomenological potentials which vary with the energy of 
the system proved useful in atomic, molecular as well as nuclear and particle 
physics. Quite often, a certain energy-dependence of the parameters results 
from their physical meaning ~ the best known illustrations are provided by 
the Klein-Gordon equation (where the potential happens to be a quadratic 
function of the energy [1]) and by the Bethe-Salpeter equations (where the 
relativistic kinematics is partially incorporated in the description of the two- 
body systems [2]). Another elementary sample of such a motivation of the work 
with an energy-dependent Hamiltonian H — H{E) may be found outlined here 
in Appendix A. 

In another group of phenomenological models, the energy-dependence in 
the "realistic" Hamiltonians H = H[E) remains artificial and serves as a mere 
ad hoc simulation of certain complicated and not too well understood effects 
within a sufficiently transparent model. In the similar pragmatic simulations, 
the fit of the experimental data is often unexpectedly successful [3] . Although 
the consequent theoretical foundation of the explicit construction of the energy- 
dependent components of the interactions may be missing, the procedure of the 
(variational or numerical) fitting itself usually offers a useful guidance towards 
an optimal choice of the form oiH = H{E) [4] . In such a situation, one can only 
regret that the merits of the simplified physics are sometimes counterbalanced 
by perceivable growth of mathematical difficulties. In this sense, the energy- 
dependent forces attract attention as a purely mathematical challenge [5]. 

Prom an abstract point of view, one of the most common theoretical sources 
of the energy-dependent potentials may be sought in a tentative projection 
of wavefunctions |^) on a certain "model" subspace of Hilbert space. After 
such a projection the unitarity of the time evolution is broken of course. The 
microscopic and energy-independent Hamiltonians H are being replaced by 
their so called effective reduced forms ifg// = Heff(E) [6]. Some of the key 
technical details of such a re-definition of the system may be found summarized 
in Appendix B. Within its overall framework, the development of a sufficiently 
flexible treatment of speciflc models seems to be an urgent task. For this reason, 
we decided to re-analyze here the general Schrodinger energy- dependent bound- 
state problem 

H{EMa) ^ E^\<f>^) (1) 

keeping in mind, first of all, the new possibilities which were opened by the 
recent growth of interest in the systems where the time evolution need not 
necessarily be purely unitary [7]. 

In section 2 we shall start our considerations by summarizing a few known 
facts about the models with the general H = H{E). We shall show there 
how the major part of the difficulties introduced by the energy-dependence of 
the Hamiltonian may be formally eliminated by the introduction of a suitable 
bi-orthogonal basis. In the new language, the enhanced flexibility of the basis 
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becomes reflected by the possibility of a modification of the scalar product 
in Hilbcrt space. In this manner, one may even try to control the physical 
contents of the theory [8]. 

In section 3 we shall pay attention to the demonstration of the applicability 
of the formalism of section 2 to a specific class of the energy dependent Hamil- 
tonians. We decided to pay attention to the next-to-solvable potentials which 
form the so called quasi-exactly solvable (QES, [9]) family. A sketchy review of 
their basic properties is presented as a separate Appendix C. The main conse- 
quence of our present energy-dependent re-interpretation of the popular QES 
Hamiltonians may be seen in the estabhshment of their new role of a source of 
non-orthogonal bases in Hilbert space. 

In the short summary of our effort in section 4 we shall emphasize that our 
Appendices A and B stressed the immediate physical and pragmatic appeal 
of the energy-dependence in if = H{E). Our subsequent choice of the QES 
illustrations of Appendix C may then play the role of their sophisticated and 
more fiexible mathematical complement, showing that although our present 
interpretation of the effective and other energy-dependent Hamiltonians seems 
to be a fairly universal recipe, its various implementations may differ in the 
degree of the insight they are able to offer. In this sense, any continuation of 
their study would be welcome. 

2 Mathematical consequences of the energy- 
dependence of H{E) 

In the spirit of the recent review [10], many (though not all of the) mathe- 
matical puzzles related to eq. (1) may be clarified when one contemplates an 
auxiliary family of the Hermitian Hamiltonians H{z) which depend on a real 
parameter z. Let this value be temporarily disentangled from the energy eigen- 
values. Then, all the Hamiltonians H = H[z) may be assigned their respective 
spectral representations, 

H{z)^Y.H^))En{z) {n{z)\. 

Only within such a broader framework, we impose the constraint En{z) — z aX, 
the very end of all the necessary constructions and for all the indices n. The 
latter constraint may happen to be satisfied at an empty or non-empty set of 
roots z = Zi{n), Z2{n), . . . , Zm(n){n). With their knowledge at our disposal (see 
an elementary example in Appendix A), one may pick up a suitable subset A 
of the "allowed" multiindices a — {n,i) & A and abbreviate Zi{n) — Eq, and 
\n{zi) = l^a). By such a construction we arrive at a formal solution of our 
energy-dependent problem (1) under very general assumptions and/or for the 
various particular choices of the parameter-dependence in H = H{z). 
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2.1 Non-orthogonal wave functions 

There exist several difficulties which arise in connection with the generalized 
Schrodinger equation (1). First of all, the standard orthogonahty relations 
between the separate bound states are lost. Even though the norm of each 

state may be fixed via a suitable re-scaling of ||0q;|| = {<Pa\'Pa)^ it is necessary 
to know and evaluate also all the off-diagonal non-vanishing overlaps 

{(t)a\(t>f)) ^ Ra,fi, a,peA. (2) 

Under certain (nontrivial) assumptions this matrix has an inverse which enters 
the decomposition of the identity projector 

I- E \^a){R-').{M- (3) 

The practical value of such a completeness of our basis (or, at least, its com- 
pleteness in the whole "relevant" Hilbert space) is significantly lowered by 
the non-diagonality of the matrix R. The sufficiently precise evaluation of 
this matrix and of its inverse R~^ may happen to become prohibitively time- 
consuming. The practical applicability of the completeness relations (3) with 
a non-diagonal R is further lowered in all the applications of the energy- 
dependent potentials which would rely upon a purely numerical solution of 
eq. (1). Within the framework of our present considerations, at least a partial 
return to semi- and non-numerical techniques will be advocated, therefore. 



2.2 Biorthogonal basis in the Hilbert space 

An inspection of the action of our family of the operators H{E) on the separate 
elements of our basis set of kets |0a) = \n{zi) with a — {n,i) E A reveals that 
the new operator K defined by its generalized spectral representation 

E {R-'),A<Pp\ (4) 

a,f3eA 

may be interpreted as equivalent to our Hamiltonian family whenever their 
action within our "allowed" space is concerned, 

K\(P^)^E^\(I)^). (5) 

By construction, the new operator K is energy-independent (i.e., no nonlin- 
earity is encountered) and non-Hermitian (this is the price we decided to pay). 
The latter property of K ^ is unexpected but its origin is clear after we 
abbreviate 

((0.1 = E (6) 

and re-interpret the above eq. (4) as a diagonal-type spectral representation 
of our auxiliary non-Hermitian quasi- Hamiltonian, 

K=Y: \ci>a)E^{{<j>^\. (7) 
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The simplicity and transparency of such a notation is based on the consequent 
use of the double-bra symbol for the left cigcn-vcctors. This convention intro- 
duces a basis which happens to be, by construction, bi-orthogonal, 

((0a|0/3) = (^a,/3, «, /3 G A . (8) 

Similarly, the abbreviated eq. (3), 

i-i:\^o){m (9) 

represents the maximally compact form of completeness relations. 



2.3 Pseudo-Hermiticity 

In a way paralleling the recent interest in certain specific non-Hermitian (for 

example, PT— symmetric) systems with real spectra [11], we may expect that 
there exists an invcrtiblc and Hcrmitian "auxiliary metric" operator 7] = r]'^ 
such that our new operators K ^ K'^ will satisfy the following quasi- or pseudo- 
Hermiticity relations 

K^'rj^rjK (10) 
[7]. After the insertion of an ansatz 

V= E l0a))t/a,/3((0/5|, (H) 

relation (10) degenerates to the mere algebraic equation 

E^Ua,p^Ua,pEp, a,PeA. (12) 

It implies that under the non-degeneracy assumption Ep ^ E^ for a ^ 13, 
all the off-diagonal elements of U must vanish while the diagonal elements 
Ua,a — da — da ^ remain variable. This defines the large set of the metrics 

^ = E \<Pa))da{{<Pa\ (13) 
aeA 

which are self-adjoint and invcrtiblc, 

77-^= E \<l>a) d-' M . (14) 

In this language each energy-dependent Hamiltonian H{E) may be better un- 
derstood via its 77— pseudo-Hermitian representation K using an arbitrary aux- 
iliary metric r) defined by eq. (13). Its role becomes clear after we introduce 
the Hermitian conjugate of eq. (5), (0o;| — (0a | E^. At the real energies 
E* — Ea it enables us to re- write this equation in the spirit of eq. (10), 

((0«|r/) = ((0«|r/) E^. 
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The non-degeneracy assumption immediately implies that 

= const {(palv- (15) 

We see that the knowledge of the (non-singular) operator r] is, up to a re- 
normalization (15), equivalent to the knowledge of the matrix R [cf. eq. (6)]. 
The main meaning of the use of rj follows from the fact that in the hypothetical 
Hilbert space where the metric would be given by the operator r), the operator 
K would become, in the hght of eq. (10), "Hermitian". 

3 Illustration: Partially solvable potentials 

The Hautot's QES model (30) of Appendix C is an example of the well- 
motivated choice of the quantum system where, by definition, a suitable selec- 
tion of the energy-dependent Hamiltonian H{E) = H(^p) = i7(o) + FW [with 
W — 1/r and F — F^'^^^\E)] is dictated by the feasibility of the evaluation of 
the overlaps (2). This means that the selected bound states i^l^e^^\r) remain 
represented by polynomials of a finite degree N > n. While we may drop the 
information about the value of the nodal count n and of the angular momen- 
tum i as redundant, we still have to keep in mind our choice of the subscript 
j which numbers all the eligible QES charges e = F/vj with j = 0,1, . . . , N. 
This allows us to compactify our notation, ipn^i = il'n,i,N,j l^^j); with the 
main purpose of using a suitable subset of the overcomplete family of these 
QES bound states as a non-orthogonal basis in the physical Hilbert space. 

3.1 QES basis and orthogonality- type relations 

Let us re-write eq. (30) as the following set of equations in the Dirac's bra-ket 
notation, 

H^o) \M, k) + W \M, k) FM,k = \M, k) Em (16) 

k^l,2,...,L{M), M = 0,1,.... 

In our particular example we shall prefer the one-charge-per-energy choice of 
L{M) = 1. The left and right eigenstates do not differ (or at least need not 
differ - see the non-Hermitian generalization of eq. (30) in ref. [13]). Thus, we 
may complement eq. (16) by its conjugate counterpart 

{N,j\ H^o) + Fnj {N,j\W = En {N,j\, (17) 

j = 1,2,..., L{N), N = 0,1,.... 

There is no a priori reason for an orthogonality between multi-indexed bra 
vectors {N,k\ = {A\ and ket vectors \N',k') = \b) . Their overlaps 
RA,b — {A\b) form a non-diagonal matrix (2) in general. We only have to 
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assume that this matrix remains invertible. Only in such a case we may employ 
eq. (3) and introduce the identity projector 

'= T. (18) 

As long as equations (16) and (17) share all their eigen-energies and eigen- 
charges, we may write down the following two alternative projected equations 

{N,j\ H^o) \M,k) = {N,j\M,k) Em - {N,j\ W \M,k) FM,k (19) 

{N,j\ H^o) \M, k) = Em {N,j\M, k) - Fmj {N,j\ W \M, k) (20) 
with {N,j) e Jbm and {M,k) e Jket- Their subtraction gives the constraint 

{FM,k - Fnj) {N,j\ W \M, k) = {Em - En) RiNjUM,k) ■ (21) 

This relation may be understood as the energy-dependence-related generaliza- 
tion of the usual orthogonality of the eigenvectors. All the matrix elements of 
Hq dropped out. For the so called Sturmian multiplets (i.e., within each sub- 
space where M — N), the left-hand-side expression must be a diagonal matrix 
with respect to its second indices j and k. For our present purposes we shall 
abbreviate {N , j\W\N ^ j) = wnj. All the other matrix elements of W may be 
treated as defined by eq. (21). Once we know all the overlap matrix R, just the 
diagonal elements w^j remain unspecified and, whenever needed, their values 
must be generated by an independent calculation. 



3.2 Perturbations and non-QES states 

Whenever we leave the safe QES domain and contemplate the more general 
bound-state problem (30) at a generic charge F ^ F^^^^\ we encounter a non- 
elementary (i.e., numerical or perturbative) diagonalization of Schrodinger 
equation 

'h^o) + FW\ \^) = E{F) I*)., 
Assuming that F ^ F^^^^^ and using eq. (18) we may insert 

and arrive at the infinite-dimensional linear algebraic equation 

^ {A\ H{F) \h) R,,c (C'l*) = £^ A e J^ra 

i.e., matrix equation 

Z{E,F)h = 0, (22) 

where 

ZA,b{E, F) = {A\ H^o) \b) - E {A\b) + F {A\ W \b) . 
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In the preparatory step we must evaluate all the input matrix elements. This 
step is usually the most time-consuming part of the algorithm. Fortunately, 
it can quite efficiently be optimized, in the present QES setting, by the use 
of all the available orthogonality-type identities. Thus, we recall eq. (20) and 
eliminate all the matrix elements of -f^(o)- This means that in eq. (22) the 
costly input information becomes reduced to the mere evaluation of the matrix 
elements of the Coulombic W{r) — 1/r, 

ZA,b{E, F) = {F- Fa) {A\ W \b) - {E - Ea) {A\h) . 

In the second step we keep M ^ N (i.e., we stay out of the Sturmian subspaces 
or diagonal blocks in the matrix Z) and postulate the absence of a random 
degeneracy of charges. This means FM,k 7^ F^j so that we are permitted 
to re-arrange the orthogonality-like relation (21) into definition representing a 
further vital reduction of the necessity of the excessive numerical integrations, 

(TV, j| W |M, k) = - En j^^^.^^^^^^^ ^ M^N . 

^M,k ~ -fN,j 

The final, maximally reduced form of our linear Schrodinger non-QES algebraic 
problem then reads 

, , En - Ek , 

WN,j riN,j + 2^ — — ni^N,j),{K,p) n.K,p - 

E - En ^ 
- ^ — 2^ n(N,j),{M,k) fiM,k , 

^ - M,k 

i = l,2,...,L(iV), N^0,1,... . 

Summarizing, any numerical or perturbative solution of this algebraic system 
requires just the knowledge of the overlaps R and of the single array of the 
special diagonal Coulombic matrix elements wnj- 



4 Summary 

Whenever we try to parallel the orthogonality proof as it works in the standard 
energy-independent cases [1] , we merely obtain a very formal relation [10] 

{<f>p\ [H{Ep) - H{E^)] 10,) = {Ep - E^) (0^10,) . (23) 

In spite of its comparative weakness (and in spite of its virtually negligible 
role in the purely theoretical considerations of section 2 above) , we succeeded 
in demonstrating that its contents remain non-empty and that its practical 
implications may be very useful. First of all, this "weak orthogonality" relation 
has been shown to imply that whenever the energies coincide. Eg = Ea-, the 
self-overlaps {4>j3\4>a) remain undetermined. Less trivial is the observation that 
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for the non-degenerate spectra with ^ Ea, formula (23) admits the non- 
vanishing of the overlaps i?^ q, and relates their values to the left-hand-side 
matrix elements. This is one of the key consequences of eq. (23), with practical 
merits which vary with the explicit form of dependence of the Hamiltonian 
operator H{E). This point of view has been more explicitly supported by our 
illustrative energy-dependent re-interpretation of the current and popular QES 
examples. 

In the broader context involving the generic energy-dependence of virtu- 
ally all the effective Hamiltonians we re-interpreted all the energy-dependent 
(i.e., in a way, non-linear) Schrodinger equations (1) as equivalent to the non- 
Hermitian but, by construction, fully linear algebraic eigenvalue problems (5). 
Such a step has been shown to faciliate our understanding of the theory as 
well as of the deeper mathematical nature of our original equation. The re- 
lated possible scalar-product re-interpretations may make it more relevant in 
the context of physics, along the lines discussed much more thoroughly in the 
review of quasi- Hermiticity [8] as well as in its more recent continuations [7, 12]. 

In the future extensions of our present note one could select several interest- 
ing directions. Firstly, whenever the deviation \ = F — F^^^^^ of the charges 
remains sufficiently small, one feels tempted to construct the spectra E = E{X) 
perturbatively, in the form of a power series in A. This could complement the 
existing QES-related perturbative studies based on different principles [14] . In 
a broader setting, our present recommendation of the use of unusual biorthog- 
onal bases might also lead to a direct new progress in the area of perturbation 
theory itself. 

As we already emphasized, the diagonalization of an operator which de- 
pends on its own eigenvalues is not, strictly speaking, a linear problem. In this 
sense the incompletely elementary QES class seems particularly suitable as an 
illustrative example at an introductory stage. In a way emphasized in ref. [10], 
all continuations of such a direction of development will be well motivated since 
even many apparently elementary energy-dependent interactions fail to admit 
a non-numerical treatment. Thus, a systematic classification of all the solvable 
cases would be of a paramount theoretical as well as practical importance. 

Last but not least, we should not forget that also the recently popular 
replacements of the Hermitian H = H[E) by their non-Hermitian descendants 
with real spectra [13] could offer another inspiration for a continuation of our 
present study. 
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Appendix A: Energy dependence with an origin 
in physics 

Under the influence of experimental data, even such a certainty as the energy- 
independence of the mass of a particle may require a critical re-evaluation. 
Such a critique finds its support in the decays of mesons where, tradition- 
ally, the vector-meson dominance offers a parameter-free explanation of the 
overall character of the data [15]. Still, the standard choice of the form factor 
F(t) of the p meson leads to a perceivable discrepancy between the parallel 
descriptions of the decays 7r+/i+/i^ and Tr+e+e". The remedy 

of this discrepancy has been found [16] in the variability of the dilepton mass 
M — tphys which is perceivably different in the above two processes. In accord 
with Isgur et al [17], one must work with the running mass — > rh^{t) in 

2 

Tfl 

Fit) = . (24) 

The variability of the running mass with the cross-channel energy t is necessary 
for a consistent interpretation of the analyticity of the propagators. Explicit 
calculations lead to the further modification of eq. (24) and an imaginary shift 
appears in 

^^^^^ ml{t)-t-impV{t) ■ ^^^^ 

The determination of the "realistic" running-mass function w?p{t) represents 
the main and challenging theoretical problem [18]. In the phenomenological 
considerations of ref. [16], a satisfactory agreement between the experiment 
and its description has been achieved via a selfconsistent determination of the 
energy dependence of the mass m^(t) near the phenomenologically relevant 
dilepton invariant energies t a; En+^- and t ~ E^+e- ■ 

For a schematic clarification of some of the basic features of the above- 
mentioned energy dependence, the most elementary explicit toy model may 
be used and studied. Once we contemplate the Schrodinger equation in one 
dimension (in units h — 2) with the harmonic oscillator interaction, 

1 d^ 

— —^(r) +r^^(r) = E^(r) (26) 

m[E) dr^ 

the effects of the variability of the mass may be mimicked by an arbitrary sim- 
ulation of its energy-dependence. For the most elementary illustrative ansatz 

m{E) ^A'iE- E,f , (27) 

a re-scaling of eq. (26) leads to the new bound-state problem with the spectrum 
determined by the closed formulae 



^ = Eo± J El - n = 0,l,... ■ 
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The latter two sets are finite and exist only for ^£^q > 4. The new levels 
emerge at each new Umax = entier[{A — 4)/8]. The message of this test is 
quite persuasive - the structure of the spectrum may be modified thoroughly 
even by a very innocent-looking energy-dependent term. 



Appendix B. An artificial energy dependence 
originating from the model-space projections 

One of the most usual approaches to the realistic Schrodinger equation is vari- 
ational. Typically, people start from a suitable microscopic Hamiltonian. Such 
a quantum model is usually constructed on the basis of the correspondence 
principle. As a rule, it is energy-independent and complicated, its handling 
proves time-consuming and its numerical diagonalization appears distressingly 
slow. These difficulties force us to reduce the Hilbert space to a smaller sub- 
space. This means that we must replace the original microscopic bound-state 
problem by its model-space version or reduction. As long as our attention is 
merely paid to the certain "most relevant" subsets of all the possible degrees 
of freedom, an abstract algebraic reformulation of the above statement may be 
based on a split of the identity operator / into the projector P (on the man- 
ageable subspace) and its complement Q — I — P. The exact and complete 
Schrodinger equation H \ ^) — E\^) acquires the two- by-two partitioned form 

P{H - E)P\m) + P{H - E)Q \m) = 0, 

Q{H-E)P\^) + QiH-E)Q |^) = 0. 

The energy-dependence emerges when we eliminate the component Q |^) of 
the wave function from the second row and insert it in the first equation. The 
resulting explicit form of the reduced problem. 



P{H -E)P\^) + PHQ 



Q 



QHP\^) = 



Q{H-E)Q_ 

may be understood as a projected or "cfi'cctivc" Schrodinger equation 

Heff\^eff)^E\^eff). (29) 

The value of the energy remains unchanged while the wave function is reduced 
to a subspace, |^e//) = In general, the energy dependence appearing 

in the effective Hamiltonian Heff{E) is very complicated. 

Appendix C. Partial solvability re-interpreted 
as an energy-dependence 

For the sake of definiteness, let us pick up the two most elementary and popular 
QES examples, viz., the Hautot's [19] shifted harmonic oscillator defined at 



12 



certain exceptional charges F only, 



V'n/(r) = En/lpnA'^) (30) 



and the sextic oscillator of Singh et al [20] with certain exceptional "spring" 
constants A, 



+ ^ — r-^ + Ar^ + ar* + r'' 



^n,^'^) = en,^0n,^(?^) • (31) 



Both these models may be made mathematically completely equivalent via a 
suitable change of the variables [21] but each of them plays a slightly different 
role in physics. We may index their bound states in the same manner, viz., by 
their angular momenta i and by the number n of nodes in the wave function. 

In both cases, the essence of the QES construction lies in the requirement 
that the Taylor series for the bound-state wave functions terminate and become 
proportional to a polynomial of degree N. In the former, Coulombic case 
(30) this implies that the exceptional QES energy becomes fixed and remains 
unique and charge-independent. Thus, the spectrum will be numbered by N 
and coincides, incidentally, with the equidistant energies of the pure harmonic 
oscillator. In contrast, at each there exist as many as N + 1 different QES- 
compatible values of the charge F — -F(Arj), with j — 0,1, . . . , N (see [19] for 
details). At each energy, the admissible value of the charge becomes unique 
only after a particular choice of the index j = Jq which may vary with the 
changes of the second index N. Thus, the QES charge may be understood as 
an energy-dependent quantity, F = Fqes{E). 

In the latter, sextic illustrative example (31), the roles of coupling and 
energy become interchanged [22]. At any angular momentum i, the main 
quantum number A^ now counts the eligible couplings A = A^. The related 
(A'-l-l)— plets of the admissible energy values £„j must be generated by the spe- 
cific algorithm of ref. [20]. Although the resulting overall energy- dependence 
pattern remains very similar, its practical aspects become less comfortable since 
up to the exceptional large— £ limit [23] , the correspondence between E and A^ 
acquires an ackward numerical character at the larger N. 
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